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A family of maximally superintegrable systems containing the Coulomb atom as a special case 
is constructed in n-dimensional Euclidean space. Two different sets of n commuting second order 
operators are found, overlapping in the Hamiltonian alone. The system is separable in several 
coordinate systems and is shown to be exactly solvable. It is solved in terms of classical orthogonal 
polynomials. The Hamiltonian and n further operators are shown to lie in the enveloping algebra 
f^*) , of a hidden affine Lie algebra. 
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I. INTRODUCTION 

We shall consider the stationary Schrodinger equation 
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HtP = E^, H = --J2 T ^ + V(x 1 ,...,x n ) (1) 

O ■ <=i 1 



o 

T— I 

o 

Oh 



X 



in an n-dimensional Euclidean space E n . By analogy with classical Hamiltonian mechanics, such a system is called 
integrable if there exist n — 1 algebraically independent linear operators X a satisfying 



[H,X a } = 0, [X a ,X b }=0, a,b = l,...,n-l. (2) 
The system is called "superintegrable" if there exist k further operators {Yi, . . . , Yfe} commuting with the Hamiltonian 

[H,Yj], j = l,...k, (3) 



such that the set {H, X\,..,, X n , Y\, . . . , Y^} is algebraically independent. Note that the additional operators Yi need 
not commute with the operators X a nor amongst each other. The number of additional operators satisfies 

l<jfc<n-l. (4) 

For k = 1 we call the system "minimally superintegrable" , for k = n — 1 it is "maximally superintegrable" . 

The best known superintegrable systems in E% (and also in E n for any n > 2) are the harmonic oscillator and the 
hydrogen atom (or Kepler system in classical mechanics). The harmonic oscillator Jji-superintegrable because of the 
su(n) algebra of first and second order operators commuting with the Hamiltoniar£ra. The hydrogen atew in E n is 
superintegrable, because of the o(n + 1) Lie algebra of linear operators commuting with the Hamiltonianu 0. In both 
cases it is possible to choose different subsets of n operators commuting with each other and overlapping only in the 
Hamiltonian. Each subset corresponds to the separation of variables in the Schrodinger equation in a different system 
of coordinatesa. 

Characteristic features of these two superintegrable systems are: 

1. In classical mechanics all finite (bounded) trajectories are periodic. Moreover, Bertrand's theoremB0 tells us 
that 7/V and jr 2 are the only spherically symmetric potentials for which all finite trajectories are periodic. 
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2. In quantum mechanics these two systems are exactly solvable: their energy levels can be calculated algebraically, 
as can the degeneracies of these levels. Their eigenfunctions are polynomials in the appropriate variables, 
multiplied by some overall factor. 

3. These systems are extremely important in physical applications, both in classical and quantum physics. 

It makes sense to search systematically for superintegrable systems in classical and quantum mechanics, specially 
for maximally superintegrable ones. It can be safely assumed that they will all have the above properties 1 and 2 and 
hoped that they will also, to some degree, share property 3. 

In searches for superintegrable systems restrictions are imposed on therJorm of the commuting operators X a and 
Yi. A systematic search in E<x and E% was conducted some time agcOila. The restriction was that all operators 
involved .should be at most of second order. All superintegrable systems satisfying this restriction in Ei and E3 were 
foundtlrta. Four classes of them exist in E2, 5 maximally superintegrable (2n —1 = 5 operators commuting with H) 
and 8 minimally superintegrable ones (n + 1 = 4 operators) in E3 . These respite have been recently extended to two 
and three dimensional spaces oLconstant curvature and to complex spacestlrEj and also to certain two dimensional 
spaces of nonconstant curvatureEl. 

With the restriction to second order operators all superintegrable systems turned out to be multiseparable, that is, 
separable, in at least two different coordinate systems. In two dimensional spaces they also turned out to be exactly 
solvableEi By this we mean that their energy spectra can be calculated algebraically (by solving algebraic equations 
only)E3 til. It was also shown that superintegrable systems are obtained by considering non-Abelian algebras of 
generalized Lie symmetriesES. 

The purpose of this article is to consider a family of integrable systems in n dimensional Euclidean space for any 
n. The family, containing the n dimensional hydrogen atom as a special case, is introduced in Section 2, together 
with a set of 2n — 1 algebraically independent operators, commuting with the Hamiltonian. In Section 3 we solve 
the Schrodinger equatian-in parabolic and spherical coordinates and show that it is exactly solvable in a precise 
and well defined senseEaEa. Finally, in Section 4 we introduce parabolic rotational coordinates in E n and solve the 
Schrodinger equation in these coordinates and also in spherical ones. We also prove the exact solvability in this case. 
Some conclusions are drawn in Section 5. 



II. A FAMILY OF MAXIMALLY SUPERINTEGRABLE SYSTEMS IN E N CONTAINING THE 

HYDROGEN ATOM 



Let us first consider the hydrogen atom in n dimensional Euclidean space E n 

H = - 1 -A-l * = r=(*? + ... + 4) 1/2 . (5) 



This Hamiltonian commutes with n(n + l)/2 linear operators, namely 

d d 
Lik = Xi- Xk——, 1 <i< k <n 

dx k oxi 



1 x 

M = - y^iPaLia + L ia p a ) + -f— l<i<Tl (6) 

2 r 

0=1 

with pi = d Xi . The operators Ln~ correspond to angular momentum, A, to the n dimensional Laplace-Runge-Lenz 
vector, characterizing the Coulomb or Kepler problemoli3. Only 2n — 1 of the operators {H, Lik, Ai\ can be, and are, 
algebraically independent. They satisfy the commutation relations 

[H,L ik ] = [H, Ai\ = 

[Lij , L a b] = Sj a Lib + SibLj a — 5i a Ljb — SjbLi a 

[Lij,A k ] = SjkA, - SikAj (7) 



The commutation relations (0) in general correspond to a Kac-Moody algebraEaO. For a fixed energy, H = E they 
correspond to the Lie algebra of the rotation group 0(n + 1), the Lorentz group 0(n, 1) and the Euclidean group 
E(n) for E < 0, &> and E — 0,rfespectively. These symmetries for n = 3 were discovered implicitly by Pauliu and 



explicitly by Focka and BargmannB. 
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According to the operator approach to the separation of variablescSTliJ, separation of variables in Schrodinger 
equation is achieved by looking for eigenfunctions of a complete set of n commuting second order operators 
{H, Xi, . . . , X„_i} 

H^ = E^ X a ^^X a ^, a = 1,.. (8) 

The operators X a will be at most linear in Ai and bilinear in L^. If more than one inequivalent set of commuting 
operators exists, the system is multiseparable, i.e., separable in more than one coordinate system. 

In view of the commutation relations (Q) any set of commuting operators {Xi} can contain at most one operator 
involving A^. 

n 

E«^° ( 9 ) 

i i,k,j,m i—1 

The complete sets of commuting operators can be classified under the action of 0(n); in particular we can rotate and 
normalize so as to have a n = 1, au — for k — 1, . . . ,n — 1. Here we just give the example of the case n = 3. It is 
easy to verify by a direct calculation that in this case, precisely four inequivalent sets exist: {H, Xi, X 2 } with 

X 1 = A 3 , X 2 = L\ 2 (10) 

X 1 = A a + a(L 2 12 + Ll 3 + Ll 1 ), X 2 = h\ 2 (11) 

Xi = L 12 + L 23 + L 31 , X 2 = L\ 2 (12) 

Xi = L 12 +L 23 +L 31 , X 2 — L% z + fL 31 (13) 

They correspond to the separation of variables in parabolic rotational coordinates, shifted spheroidal coordinates, 
spherical coordinates and spheroconical coordinates, respectively. 

In each coordinate system it is possible to add further terms to the potential — 7/r in such a manner that the 
Schrodinger equation still separates. The system will remain integrable and the corresponding operators X\ and X 2 
will only be modified by the addition of a scalar function, It is also possible to preserve superintegrability and to 
require that the extended potentials should allow separation of variables in at least two coordinate systems. 

Here we will be interested in the most general potential allowing separation of variables in the same four coordinate 
systems, as the hydrogen atom itself. In E 3 there is, up to equivalence, only one such Hamiltonian, namely (see 
Ref-laW) 



1 A 1 0i 0: 



^ = -o A -7 + 3 + 5 ( 14 ) 



One triplet of commuting operators for the Hamiltonian (|1J) consists of 

1 < r , r , - r , r N , n ( 7 ft # 



X = -(p 1 L 3 i + L31P1 +p 2 L 32 + L 32 p 2 ) + 2x 3 



2. 

2r x\ x\ 



Z = L? 2 -2r 2 f§ + ^) (1 r„ 



'-I 



Another triplet can be chosen to be H and 



2 0.2 (0i , ft 



Y\ — L\ 2 + + L 31 — 2r 

Y 2 = Lj 3 ~ 20 2 ^±^ (16) 
x 2 

It is the set of 5 algebraically independent operators {H, Xi, X 2 ,Yi,Y 2 } which guarantees that the Hamiltonian ( |l4| ) 
is maximally superintegrable. 

The generalization to the n dimensional Euclidean space E n is immediate. Thus, the Hamiltonian will be 

1 a 
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with A and r as in Equation (||) . One of the two different complete sets of commuting operators can be chosen to be 
H and 



1 n ~ 1 / " _1 /3- 
X = - ^2(L nk p k + p k L„k) + 2x„ I ^ "T 

k=l \ T i=l Xi , 



*«= E ^- 2 fE^fE§V l<i<«-2 (18) 



l<i</c<i + l \i=l / \fc=l k , 



Another complete set of commuting operators is again H and 



y P = E ^- 2 (E* 2 I (E§l- i<p<»-i ( 19 ) 



p<i<k<n \^—P J \k—p K 



The two sets ( fl8| ) and (|L9|) are disjoint. If we set /3j = 0, 1 < i < n — 1, then the operator Zi will be a Casimir 
operator of the group 0(7 + 1) acting on the coordinates {xi, . . . , xi+i}. The operator Y p will be a Casimir operator 
of 0(n + 1 — p) acting on the coordinates {x p , . . . , x n } 

It is the Hamiltonian ( |l7| ) that we shall study in the following sections, first for n = 3, then for arbitrary n. 



III. EXACT SOLVABILITY OF THE SUPERINTEGRABLE SYSTEM FOR N = 3 



A. Solution by separation of variables 

Let us first consider the Hamiltonian ( |i"4| ) and the complete set of commuting operators ([l5|). We arc looking for 
eigenvalues and common eigenfunctions of the systems: 

Hip = Eil), Xif> = Xtp, Zip = kip (20) 

To do this we introduce parabolic rotational coordinates, putting 

x\ = p,ucos<f>, X2 = iAVsm<j>, x^ = 7j(^ 2 — ^ 2 ) ( 2 ^) 



In these coordinates the operators in (|20|) are 



1 / d 2 1 d d 2 1 d „ 

H = 1 1 1 h 47 

2(/i 2 + v 2 ) \dfi 2 fid[i dv 2 v dv 

1 / 8 2 2/3! 2/? 2 



2[i 2 v 2 \d(f) 2 cos 2 (f> sin 2 1 

1 / 2 ( d 2 1 d\ 2 [ d 2 li) 



/i 2 - 1/ 2 / 8 2 20i 2f3 2 



Z = 



2y?v 2 \d<t> 2 cos 2 4> sin 2 . 
d 2 20i 2f3 2 
dcf) 2 cos 2 ^ ^ ; " 2 



sin 

We see immediately that the variables separate and we can solve the corresponding ODE's to obtain 



ft = \vi(jPi - 1)) m = 2J + ;p 1 +j?2 



(22) 



i>N lt N 2> j = (sin^(cosff 1 H m e^ & ' 2+ '' a) x 

p ( P2 -i/ 2 , Pl -i/2) (cos 20)i „ (^Z^^ 2 )^™ (V^B^ 2 ) (23) 

where P l f^\z) and L™{x) are Jacobi and Laguerre polynomials respectively. We have put 
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and the eigenvalues in Equation ([20|) are equal to 

E 
X 



2{N 1 +N 2 + 2J + Pl + P2 + l) 2 
7(^1 - N 2 ) 



Ni + N 2 + 2 J + Pl + p 2 + 1 
k = -m 2 = -(2J + pi +P2) 2 



(24) 



We see that the bound state energy is given by a shifted Balmer formula and the only effect of the /3j /x 2 terms in 
the potential is to add a constant p\ +P2 to the principal quantum number. The solutions (23) are square integrable 
and correspond to bound states when J, iVi and N2 are integers. They are polynomials multiplied by a factor which, 
however, is not "universal" . It depends on the energy E and also on the angular quantum number J (since we have 
tn = 2J -\- p\ + p 2 ). 

The second set of commuting operators, namely ( |l6| ) also corresponds to the separation of variables, this time in 
spherical coordinates, chosen as 



xi = r cos f 



X2 = r sin 9 cos a, X3 — r sin u sin a 



(25) 



In these coordinates we have 

Htjj 

Y 2 ^ 
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" d 2 


2 d 
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dr 2 


r or 


Hi. 
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v cos 2 9 



8P_ 

de 2 
d 2 

da 2 



COtd d9- 



2/3: 





d 2 


r 2 ( 






w +cot 




02 y 


sin 2 


9 cos 2 a J 


2/3i 


-+ ^ ' 


cos 2 


9 sin 2 9 



^89 



1 



sin 

ip = Eip, 



' da 2 



(26) 



TT^7 5— ]1p = k 2 1p 



The coordinates separate and we obtain ip — R(r)F(9)G{a) where F and G can be expressed in terms of Jacobi 
polynomials and R in terms of Laguerre ones. 

The explicit expression for the eigenfunctions in this type of spherical coordinates is: 



*I>N, .h,h (r, 9, a) = r mi - 1/2 e- ^ 2Er L 2 ^ 11 (2V^2E r)(sin 6) mi (cos 9) Pl (cos a) P2 x 



p jm 2lPl -i/2) (cog 20)p j-i/2, P2 -i/2) (cqs 2a) 



with eigenvalues equal to: 
E = - 



7 



2(AT + 2 Ji + 2 J 2 + pi + p 2 + 1) 



(27) 



(28) 



fc 2 = -? 



mi = 2Ji + 2J 2 +pi+p 2 + -, 



m 2 = 2J 2 + p 2 



B. Exact solvability and underlying affine Lie algebra 

We have established that the potential (Q) provides a Hamiltonian that is maximally superintegrable and mul- 
tiseparable. Let us now turn to the question of exact solvability. A Hamiltonian is exactly solvable if its spectrum 
can be calculated algebraically. This occurs if it can be explicitly transformed into block diagonal form where each 
block is finite dimensional. This in turn means that there exists an infinite flag-of_finite dimensional subspaces in the 
Hilbert space Ti of bound state solutions that is preserved by the HamiltonianETo 

H 1 cn 2 c---cn, HH^H^ (29) 

Typically this occurs under the following circumstances 
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1. The bound state wave functions are polynomials in some variables, possible multiplied by some factor <?, i.e., 
-0 = gP. We then have 

h = g~ x Hg, hP = EP, (30) 

that is, there exists a gauge transformation and a change of variables to a new Hamiltonian h that has polynomial 
eigenfunctions. 

2. The gauge transformed Hamiltonian h is an element of the enveloping algebra of some affine Lie algebra L, a 
basis of which can be realized by the operators 

d d 

Ki = —, M ik = Si —, i,ke{l,...,n} (31) 
osi ds k 

(in some coordinates Si). 

Let us now investigate the Hamiltonian (14) and the commuting set of operators (]l5|) and ( |l6| ) from this point of 
view. 

Consider the wave function (E3T) in parabolic coordinates. They do have the required form 



i>Ni,N 2 ,J = 9PNi,N 2 ,j(s,t, z), s = V-2Efi 2 7 t — \J~2Ev 2 , z = cos2</> (32) 

-(s+t)/2 



i + z\ pi/2 /i-^y 2/ V st w 2 



where P is a polynomial (a product of three polynomials in one variable each). To proceed further we must get rid 
of the conformal factor (p 2 + v 2 )^ 1 figuring in equation (22) and replace the system ( p0| ) by 

Qo*p = 2 7 V>, Qiip = 2Xi/j, Zip = -m 2 V> (33) 



Qo = (/i 2 + v 2 ){H -E) + 2j 

1 / d 2 Id d 2 Id 

2 \ d/i 2 [i dfi dv 2 v dv 
Qi = 2X + (/i 2 - v 2 ){H-E) = 

1 / d 2 Id d 2 Id 

2 \ d[j, 2 [i d\i dv 2 v dv 

Q 2 = Z. 



\ m 2 






) + - 


(m 2_ 





\ m 2 






) + - 


L 2 


v 2 ) 



E(p 2 - v 



2\ 



We see here a phenomenon which has been called "metamorphosis"EJ or "migration"E3 of the coupling constant. In 
equation (^) the energy E plays the role of the frequency of a harmonic oscillator whereas the Coulomb coupling 
constant 7 plays the role of an eigenvalue of Qq. The other eigenvalues, A and m 2 , remain eigenvalues (of Qi and Z 
respectively) . 

Similarly as in the case of potentials containing the Coulomb atom as a special case in 2 dimensions, it is the system 
( |33| ) (rather than the original system fl20|)) that is exactly solvable in the sense defined above. Indeed, let us gauge 
rotate the operators Qo, Qi and Q2 and transform to the variables s, t and z. We obtain 



Qo = g^Qog = -2V^2E (sd 2 + (m + 1 - s)d s + td 2 + (m + 1 - t)d t - m - l) , 
Qi = g^Qig = -2V^2E (sd 2 + (m + 1 - s)d s - td 2 - (m + 1 - t)d t ) , (34) 
<? 2 = g~ x Qig = 4(1 - z 2 )dz 2 + 4( Pl -p 2 - ( Pl + P2 + l)z)d z - ( Pl + P2 ) 



2 



We see that Q^, fi = 0,1,2 lie in the enveloping algebra of the direct sum of three affine Lie algebras, saff (1,11) 
50^(1,^ © saff (1,11), realized by 

{d s ,sd s ,dt,td t ,d z ,zd z }. (35) 



Now let us consider the two remaining operators Yy and Y% of (j2^). They are not diagonal in the basis that we use 
(where iJ, X, Z and equivalently Qo, Qi, Q2 are diagonal). They do however commute with the Hamiltonian so they 
can only mix states of equal energy. We also have 



G 



[Y 1 ,Z}=0, [Y 2 ,Z]^Q, [Y 1 ,X}^0. (36) 

It follows that Y\ also preserves the quantum number k = —m 2 of equation (|20|), whereas Y 2 mixes all states of 
a given energy. The gauge factor g depends not only on the energy, but also on m, the eigenvalue of Z. Thus, 
Y\ = g Y\g should transform polynomials into polynomials, whereas Y2 = g~ x Yig is not obliged to. Performing the 
gauge transformation and change of variables we find 

Yi = g~ x Y x g = st{d s - t ) 2 - (m + l)(s - t)(d s - d t ) - m{m + 1). (37) 

The algebra underlying this expression includes td s and sdt, in addition to the elements listed in ([35|). We recognize 
this to be the Lie algebra saff (2, R) © saff (1, R) 

Superintegrable systems, including the hydrogen atom as special case, in 2 dimensions were associated with the 
algebra saff (2, R). The extension to n = 3 is seen to lead to saff (2, R) saff (I, II), not to saff (3,R) as one might 
have expected. It follows from expression (|37]) that Y\ will take polynomials into polynomials and indeed we have 

%P Nu n 2 ,j = -[2N X N 2 + (JVi + N 2 + m)(m + 1)}P Ni ,n 2 ,j + 

(JVi + m)(N 2 + l)P Nl _ ltNa+liJ + (JVi + 1)(JV 2 + m)P Wl +i,JV 2 -i,j. (38) 

The operator Y2 does not take polynomials into polynomials and cannot be written as an element of the enveloping 
algebra of an affine Lie algebra. 

The one-dimensional equations appearing above are easily show p i tQ ibe related to the standard types in the classi- 
fication of exact and quasi-exact solvable one-dimensional systemsr^r 6 !. 



IV. EXACT SOLVABILITY OF THE SUPERINTEGRABLE SYSTEM IN E N 

A. Solution by separation of variables 

We now consider the Hamiltonian ( [T^ ) for arbitrary n. It allows separation of variables in many coordinate 
systems. We shall use parabolic rotational coordinates corresponding to the set of operators ( [Lsl ) and spherical ones 
corresponding to the set (|l9|). The parabolic coordinates (/x, v,Q\,... Qn^-i) are defined by the relations 

X\ = [IV COS 81 COS 6*2 . . . cos6*„_ 3 cos6>„_ 2 , 
x 2 = [iv cos B\ cos 82 ■ ■ ■ cos # n _3 sin 9 n -2, 
Xs = l±v cos 6*i cos 62 ■ ■ ■ sin 6* n -3 , 



(39) 



X n -2 = [iv costq Sint>2, 
x n-\ — fivsmdi, 

x„ = i(/i 2 - v 2 ). 

We put /3i — Pi(pi — l)/2 in equation ([l7|). The eigenvalue problem that we have to solve is 

Hip = Etp, Xi/j = Xij}, Zii> = kii>, 1 < / < n - 2. 



(40) 



In parabolic coordinates the operators H and X will involve all variables, the operators Zi will involve the angles 
only. Indeed, we have 



H = - 



1 



d 2 



-2d d 2 



-2 d 



2(fj? + v 2 ) \dj-i 2 [i dfi dv 2 v dv 
Pi(Pi - 1) P2(P2 - 1) 



2l 



1 



\i 2 + v 2 2[f?-v 2 

Pn-lijPn-l - 1 



n-2) 



cos 2 61 ■■ ■ cos 2 6 n -2 cos 2 61 ■ ■ ■ cos 2 #„_ 3 sin 9 n - 2 



sin 2 6»i 



(41) 



X 



1 



2{n 2 + v 2 ) 



d 2 



n-2 d 



^d[i 2 fi dfi J ^ \dv 2 



d 2 



n-2 d 
v dv 



A(5„_ 2 ) 



22 22 
[i A — v z M v 

^ f-l 2 + V 2 2jJL 2 V 2 

Pljpl - 1) _ 

cos 2 &i ■ ■ ■ cos 2 9 n - 2 cos 2 9i ■ ■ ■ cos 2 6»„_ 3 sin 2 



P2{P2 - 1) 



Pn-l{Pn~l - 1) 



sin 2 6»i 



(42) 
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where A(S' n _2) is the Laplace operator on an n — 2 dimensional sphere 
The operators Zi satisfy 



Z ^ = ^2 T~n i> = hip, 



if 

del 



d 2 


PiiPi - 1) 


P2{P2 


-1) 


901-2 


COS 2 9 n - 2 


sin 2 t 


»»-2 


d 


+ kl 


P3{P3 


-1) 




COS 2 6»„_ 3 


sin 2 t 


'n-3 



= I tan6»„_ 3 — h 9 r - 2 -^ ] r/> = fe^i 



'n— 3 

and in general 

^ = fa/ - l)tantf w _,_ lM ^— + "I' 1 - P '^ +1 ' 1) ) * = W (43) 

V^-j-i w B _n cos 2 6»„_ ( _i sin^^i-i / 

l<Z<n-2, fco = -piCpx - 1). 

We write 

n-2 

V> = M{ii)N{v) FJ FiiQn-i-t) (44) 

and solve (|43"|) to obtain Fx in terms of Jacobi polynomials 

Fi(*„-J-i) - (sin0 n _ i _ 1 )w+ 1 (cos0„_ / _ 1 ) m - 1 + 1 -'/ 2 pjf +1 - 1/2 ' mi - l) (cos20„_ ; _ 1 ) (45) 



//^ = 2^J i +^p i + ^, /cz-L^-m 2 (46) 

= 1 8=1 



The equations for M(fi) and iV(i/) are obtained from equation (41) and (|42j) once the angular part is replaced by 
k n -2- The final result is that the wave functions are: 



n-2 

II ( sin 0n-i-i) Pl+1 (cos 9, M ) m ^ +1 - l / 2 L^- 2 {V^2E^ 2 )L^;- 2 (^lEv 2 ) x 



n-2 



n^r^'^W^-i), <r = 2^J i + 5> (47) 



,(p l+1 — 1/2, n»i_i) 

Ji 

i=l i=l i=l 

The energy is given by a shifted Balmcr formula 

7 2 

E = 4- i (48) 

2(^ + ^ + 2^=1 Ji + YZi ft + ^f 1 ) 2 

and the remaining quantum number is 

A = 1^1 ~ ^ (49) 

N ± + N 2 + 2 YZZi Ji + TZi Pi + ^ 

We see that the case of n arbitrary is a straightforward generalization of n = 3 and involves the same functions, 
namely, Jacobi and Laguerre polynomials. 

Obviously, one can also solve in spherical coordinates. In fact, formulas (EM) can be written as 



I n— 1 

x a =[ivs a , x n = -(^ 2 - v 2 ) o = l, ...,n-l, ^\s 2 = 1 (50) 



a=l 



and we could introduce any coordinates on the S n -2 sphep©ij&atallow separation of variables in the Laplace-Beltrami 
equation. For a discussion of such coordinate systems seeEjEltSl. 
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We will write for the sake of completeness the explicit expression of the eigenfunctions in the following set of 
spherical coordinates on the S n -i sphere (which are a generalization to dimension n of those we used in the case 
n = 3, see Equation (EE 



x\ = rcosc/i 

X\ = r sin 0i cos #2 

x n -i = r sin 0i • • • sin0„_2 cos0„_i 
x n = r sin 9\ ■ ■ ■ sin 9 n -i sin n -i 



and the Hamiltonian can be written in these coordinates as 



H 



r r 



1 

2T 2 



dl + (n-2)cot9 1 de 1 - Pl{P \ Q 1) 

1 COS 2 0\ 



sin 2 0i 



d 6 2 + (n - 3) cot 2 <9e 2 



P2(P2 - 1) 
COS 2 02 



sin t/„_ 3 



1 



d 8„ , + COt0„_ 2 9fl _ 2 



P»- 2 (p»-2 - 1) 
COS 2 0„_2 



+ 



COS 2 0„_i 



The set of Y;, Z = 1, . . . ,n. — 1 operators are: 

y ; = a 2 ; + (n-z-i)cot0id ei 

v c,2 Pn-lOn-1 - 1) 

and the eigenvalue equations: 



pi{pi - 1) fc. 



-5 , {= 1, ... 71 — 2 

COS 0; sin 0; 



COS 2 n _i 



HtP = EtP, Y l G l {6 l ) = k l G l (e l ), Z = l,...,n-1, i\> = R(r) JJ G,(0,) 

1=1 

can be easily solved. The solution for the angular part is (m„ = —1/2): 

Tl—l n— 1 

n = n( sin ^) m!+i+i_(n ~' )/2 ( cos ^) wp jr +i,pi ~ 1/2) ( cos2 ^) 



z=i 



and for the radial part: 
The energy is written as: 



R(r) = r 



— _mi-(n-2)/2 -V=2£r r 2mi 



2(AT r + 2 Elli Ji + ElLi Pn-i + |(n - l)) 2 
and the eigenvalues of the operators Yi are: 

^ n—l n—1 _^ 

fc; = -(n-Z-l) 2 -m 2 , mi = 2^ Jj + + -(n - / - 1), J = l,...,n-1 

Finally, the eigenfunctions are: 

^,j 1> ...,j n _ 1 (r, 0i, • ■ • , 0„-i) = ^-(-2)/^-^^^ (2 V^2£r) x 



n— 1 



n 



(sin0O m,+1+1_( "" O/2 (cos0O J 

1=1 L 

p(m !+ i,p ( -l/2) ( 



'(COS 20;) 



(51) 



(52) 



(53) 
(54) 

(55) 

(56) 
(57) 

(58) 



(59) 
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B. Exact solvability 



The exact solvability of the system ( p7| ) for general n can be treated in the same way as for n — 3. We can gauge 
transform each of the operators in the set (|l8| ) separately and transform to the variables 

s = V^2~E/j 2 , t = v^2Ev 2 , z n - l+ x = cos20 w _,+i, l = l,...,n-2 (60) 

Before doing this, we again introduce Qq and Qi as in equation ([33]). 
The final result is 

Qo + Qi = g^iQo + Qi)g = -2y/-2E(2sd 2 + 2(1 + m„_ 2 - s)d s - m„_ 2 - 1) 

Qo - Qi = 9~ l {Qo - Qi)9 = -2^/=2£(2^^ t 2 + 2(1 + m n _ 2 - i)d t - m n _2 - 1) (61) 

Z, = g- x Z l9 = 4(1 - 4-i-i)-e4 + 

UZ n-l-l 

4 (mi-i - Pi+i + 7;- \ Pi+i + m i-i + I J z n _z_i ] d Zn _ l _ 1 + 



2 

— ^ (rnj_i +p;+i)(m ; _i + 1) 

We see that the entire set of operators {Qo, Qi, ■ ■ ■ , lies in the enveloping algebra of direct product of n 

special affine Lie algebras saff(l, R). 

Finally let us turn to the other complete set of commuting operators (|l9|), associated with the separation of variables 
in spherical coordinates. Among these operators there is just one, namely Y\, that commutes with all the operators 
Zi. We have 

[Y u Zi}=0, [Y 1 ,X]^0, [Y p ,Z t }^0, l</<n-2, 2<p<n-l (62) 
Thus, Y\ will take polynomials into polynomials but {Y2, • ■ ■ , ^n-i} w iH not. We have 

Y x = g-^Yyg = st(d s - d t ) 2 - (m„_ 2 + l)(s - t)(d s - d t ) - 

. . (n — 3)(n — ll .„„. 
m„_ 2 (m„_ 2 + l) + ^ 'j '- (63) 

Finally we see that the "hidden Lie algebra" that is not a symmetry algebra of the problem, but underlies its exact 
solvability is saff(2, R) © [saff(l, R)]i © • • • © [saff (1, R)]„- 2 generated by 

{<9 S , d t ,sd 8 , td t , sd t , td s ,d Zl ,zid Zl ,. . . , d Zn ^ 2 , z n ^ 2 d Zn _ 2 } (64) 



V. CONCLUSIONS 



Supcrintegrability and exact solvability were defined in completely different ways, though both have a group theo- 
retical underpinning. Superintegrability fot-a Hamiltonian system is defined by the requirement that there be more 
integrals of motion than degrees of freedomtll. It can be characterized by the fact that the corresponding Schrodineer 
equation allows a nonabelian algebra of generalized symmetries, containing an n-dimensional Abelian siaJpalgebraEj. 
Exact solvability is defined by the requirement that the energy spectrum can be calculated algebraicallyE3~EJ. It can 
be characterized by the fact that the -Hamiltonian lies in the enveloping algebra of a certain type of finite dimensional 
affine Lie algebra. It was conjecturedE3 that all maximally superintegrable systems are exactly solvable. In this article 
we have confirmed the conjecture for the considered integrable system in E n . 

The exact connection between superintegrability and exact solvability remains an open problem. 



ACKNOWLEDGEMENTS 



The research reported in this article was performed while P.W. was visiting the Departamento de Ffsica Teorica 
II of Universidad Complutense of Madrid. The final version was written while he was visiting the Isaac Newton 



10 



Institute in Cambridge. He thanks both Institutions for their hospitality and support. The research of P.W. was 
partially supported by research grants from NSERC of Canada and FCAR du Quebec. The research of M.A.R was 
partially supported by DGES Grant PB98-0821 of Spain. 



1 J. Jauch and E. Hill, On the problem of degeneracy in quantum mechanics. Phys. Rev. 51, 641-645 (1940). 

2 M. Moshinsky and Yu. F. Smirnov, The Harmonic Oscillator in Modern Physics. Harwood, Amsterdam, (1996). 

3 W. Pauli, Uber das Wasserstoffspektrum von Standpunk der neuen Quanten-mechanik. Zeits. Physik 36, 336-363 (1926). 

4 V. A. Fock, Zur Theorie des Wasserstoff atoms, Zeits. Physik 98, 145-154 (1935). 

5 V. Bargmann, Zur Theorie des Wasserstoffatoms, Zeits. Physik 99, 576-582 (1936). 

6 M. J. Englefield, Group Theory and the Coulomb Problem. Wiley, New York (1972). 

7 H. V. Macintosh, Symmetry and Degeneracy, in Group Theory and its Applications, edited by E.M. Loebl, (Academic Press, 
New York, 1971), Vol. II, pp. 75-144, 

8 E. G. Kalnins, W. Miller Jr, and P. Winternitz, The group 0(4), separation of variables and the hydrogen atom. SIAM J. 
Appl. Math. 30, 630-664 (1976). 

9 J. Bertrand, Theoreme relatif au mouvement d'un point attire vers un centre fixe. Comptes Rendus Ac. Sci 77, 849-853 
(1873). 

10 H. Goldstein, Classical Mechanics (Addison- Wesley, Reading, MA, 1990). 

11 I. Fris, V. Mandrosov, J. Smorodinsky, M. Uhlff, and P. Winternitz, On higher symmetries in quantum mechanics. Phys. 
Lett. 16, 354-356 (1965). 

12 P. Winternitz, J. Smorodinsky, and M. Uhlff , Symmetry groups in classical and quantum mechanics. Yad. Fiz. 4, 625 (1966), 
Sov. J. Nucl. Phys. 4, 1326 (1967). 

13 A. Makarov, J. Smorodinsky, Kh. Valiev, and P. Winternitz, A systematic search for nonrelativistic systems with dynamical 
symmetries. Nuovo Cimento A 52, 1061-1084 (1967). 

14 N. W. Evans, Superintegrability in classical mechanics. Phys. Rev. A 41, 5666-5676 (1990). 

15 N. W. Evans, Group theory of the Smorodinsky- Winternitz system. J. Math. Phys. 32, 3369-3375 (1991). 

16 N. W. Evans, Superintegrability of the Winternitz system. Phys. Lett. A 147, 483-486 (1990). 

17 E. G. Kalnins, W. Miller Jr, and G. S. Pogosyan, Superintegrability and associated polynomial solutions: Euclidean space 
and the sphere in two dimensions. J. Math. Phys. 37, 6439-6467 (1996) 

18 E. G. Kalnins, W. Miller Jr, and G. S. Pogosyan, Completeness of multiseparability in E 2 ,c- J- Phys. A 33, 4105-4120 
(2000). 

19 E. G. Kalnins, W. Miller Jr, and G. S. Pogosyan, Completeness of multiseparability on the complex 2-sphere. J. Phys. A 33, 
6791-6806 (2000). 

20 E. G. Kalnins, W. Miller Jr, and G. S. Pogosyan, Completeness of multiseparability in two-dimensional constant curvature 
spaces. J. Phys. A 34, 4705-4720 (2001). 

21 E. G. Kalnins, J. M. Kr ess, and P. Winter nitz, Superintegrability in a two-dimensional space of non-constant curvature. 
(LANL preprint archives math-ph/0108015). J. Math. Phys. (to appear). 

22 P. Tempesta, A. V. Turbiner, and P. Winternitz, Exact solvability of superintegrable systems. J. Math. Phys. 42, 4248-4257 
(2001) 

23 A. V. Turbiner, Lie algebras and linear operators with invariant subspaces, in Lie Algebras, Cohomologies and New Findings 
in Quantum Mechanics, edited by N. Kamran and P.J. Olver, (AMS, Providence, 1994), Vol. 160, pp. 263-310. 

24 A. V. Turbiner, Quasi-exactly solvable problems and sl(2) algebra. Commun. Math. Phys. 118, 467-474 (1988). 

25 M. B. Sheftel, P. Tempesta, and P. Winternitz, Superintegrable systems in quantum mechanics and classical Lie theory. J. 
Math. Phys 42, 659-673 (2001). 

26 J. Daboul, P. Slodowy, and C. Daboul, The hydrogen algebra as a centerless twisted Kac-Moody algebra. Phys. Lett. B 
317, 321-328 (1993). 

27 C. Daboul and J. Daboul, From hydrogen atom to generalized Dynkin diagrams. Phys. Lett. B 425 135-144 (1998). 

28 P. Winternitz and I. Fris, Invariant expansions of relativistic amplitudes and subgroups of the proper Lorentz group. Sov. J. 
Nucl. Phys 1, 636-643 (1965). 

29 P. Winternitz, I. Lukac, and Ya. Smorodinskii, Quantum numbers in the little groups of the Poincare group. Sov. J. Nucl. 
Phys 7, 139-145 (1968). 

30 J. Patera and P. Winternitz, A new basis for representations of the rotation group. Lame and Heun polynomials. J. Math. 
Phys. 14, 1130-1139 (1973). 

31 W. Miller Jr, J. Patera, and P. Winternitz, Subgroups of Lie groups and separation of variables. J. Math. Phys. 22, 251-260 
(1991). 

32 W. Miller Jr, Symmetry and the Separation of Variables (Addison- Wesley, Reading, MA, 1997). 



11 



E. G. Kalnins, Separation of Variables for Riemannian Spaces of Constant Curvature (Longman, Burnt Mill, 1986). 

J. Hietarinta, B. Grammaticos, B. Dorizzi, and A. Ramani, Coupling-constant metamorphosis and duality between integrable 

Hamiltonian systems. Phys. Rev. Lett. 53, 1707-1710 (1984). 

A. Gonzalez-Lopez, N. Kamran, and P.J. Olver, Normalizability of One-dimensional Quasi-Exactly Solvable Schrodinger 
Operators. Commun. Math. Phys. 153, 117-146 (1993). 

M. Shifman, New findings in quantum mechanics (Partial algebraization of the spectral problem). Int. J. Mod. Phys. A 4, 
2897-2852 (1989). 

N. Ya. Vilenkin, Special Functions and the Theory of Group Representations (AMS, Providence, R.I., 1968). 

N. Ya. Vilenkin, G. I. Kuznetsov, and Ya. A. Smorodinskii, Eigenfunctions of the Laplace operator realizing representations 

of the groups f/(2), 517(2), 50(3), £7(3) and 5(7(3) and the symbolic method. Sov. J. Nucl. Phys. 2 645-655 (1965). 

A. A. Izmest'ev, G. S. Pogosyan, A. N. Sissakian, and P. Winternitz, Contractions of Lie algebras and the separation of 

variables. The n-dimensional sphere. J. Math. Phys. 40, 1549-1573 (1999). 

G. S. Pogosyan and P. Winternitz, Separation of variables and subgroup bases on n-dimensional hyperboloids (to be pub- 
lished) . 



12 



